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O To determine whether particular sources of gravitational radiation will be de- 
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tectable by a specific gravitational wave detector, it is necessary to know the 

bp. 

sensitivity limits of the instrument. These instrumental sensitivities are often 
depicted (after averaging over source position and polarization) by graphing 
the minimal values of the gravitational wave amplitude detectable by the in- 
strument versus the frequency of the gravitational wave. This paper describes 
in detail how to compute such a sensitivity curve given a set of specifications 
for a spaceborne laser interferometer gravitational wave observatory. Minor 
errors in the prior literature are corrected, and the first (mostly) analytic 
calculation of the gravitational wave transfer function is presented. Example 
sensitivity curve calculations are presented for the proposed LISA interfer- 
ometer. We find that previous treatments of LISA have underestimated its 
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sensitivity by a factor of \/3- 
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I. INTRODUCTION 



Advances in modern technology have ushered in an era of large laser interferometers 
designed to be used in the detection of gravitational radiation, both on the ground and 
in space. Such projects include the LIGO and VIRGO |],|2| ground-based interferometers, 
and the proposed LISA and OMEGA BM space-based interferometers. As these detectors 
come on-line, a new branch of astronomy will be created and a radically new view of the 
Universe is expected to be revealed. With the era of gravitational wave astronomy on the 
horizon, much effort has been devoted to the problem of categorizing sources of gravitational 
radiation, and extensive studies are underway to determine what sources will be visible to 
the various detectors. 

Typically, the sensitivity of detectors to sources of gravitational radiation has been il- 
lustrated using graphs which compare source strengths (dimensionless strain) to instrument 
noise as functions of the gravitational wave frequency. Many different types of plots have 
appeared in the literature, ranging from single plots of spectral density to separate ampli- 
tude plots for each class of source. When considering the possibility of observing a new 
source, or comparing aspects of various proposed gravitational wave observatories, it is im- 
portant to be able to generate consistent and accurate noise curves for a given instrument 
and to understand what assumptions have gone into generating the curves. This is especially 
important in the case of spaceborne gravitational observatories which are sensitive to Low 
Frequency (LF) gravitational waves (in the band from 10 -4 Hz to 1 Hz). In this LF band, 
the sources of radiation include both known continuous sources at well-determined strengths 
and frequencies |J, for which a signal-to- noise ratio only slightly greater than one is needed 
for detection, as well as speculative short-lived "burst" sources requiring a much greater 
signal-to-noise ratio for detection. 

This paper reviews the mathematical formalism and methodology for generating noise 
curves for a class of spaceborne gravitational wave interferometers. A synthesis is provided 
here of material which has hitherto been scattered across the literature, and a variety of 
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new results are incorporated. The gravitational wave transfer function of an interferometer, 
averaged over source direction and polarization, is calculated here for the first time,[] and 
a number of minor errors in the existing literature are corrected. In particular, we find 
that several previous estimates of the sky- and polarization-averaged sensitivity of space 
interferometers have understated their sensitivity by a factor of Although the results 
of this paper are applicable to any spaceborne laser interferometer system designed for 
gravitational wave detection, specifics of the proposed LISA mission are used as an example. 

The outline of the paper is as follows. In Section |l| the concept of instrumental operation 
for a spaceborne interferometer is described. Section |T| discusses noise sources that limit 



the sensitivity of the detector and Section [TV] presents the mathematical formalism for 
generating sensitivity curves from these noise sources. Section [V] considers the various 
types of sensitivity curves which currently exist in the literature and reconciles the different 
methods. 



II. INSTRUMENT OPERATION 

A. Instrument Design 

The common design concept for proposed spaceborne interferometers consists of a con- 
stellation of probes arranged in an equilateral triangle inscribed on the circle of the probes' 
relative orbits. In the simplest configuration, 3 probes are used to form a single Michelson 
interferometer, as shown in Figure [l]. The probe at the vertex of the angle corresponds to 



1 While this paper was in preparation, we became of aware of recent work by Armstrong, Es- 
tabrook, and Tinto || which numerically estimates the gravitational wave transfer function using 
Monte Carlo simulations to conduct the averaging. A visual inspection of the sensitivity curve 
derived from the transfer function derived by M shows agreement with the results derived in this 
paper. 
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the central mirror in the interferometer, while the two probes on the ends of the arms cor- 
respond to the end mirrors. The effect of a gravitational wave passing by the detector is to 
stretch or contract space in the arms of the interferometer. The effects in the two arms will 
in general be different, due to the different orientation of the two arms in space. The wave 
is detected by monitoring the times of flight of laser signals between the probes in order 
to observe and measure this difference. The basic operation of a space gravitational wave 
detector consists in measuring the phases of the incoming laser signals relative to those of 
the outgoing signals in both arms of the interferometer and differencing these relative phases 
to cancel common phase noise in the two arms. 

In order to avoid spurious motions of the spacecraft that could mimic the effect of a 
gravitational wave on the laser tracking signals, each spacecraft is equipped with a position 
control system. This system consists of, first, an accelerometer that measures the motion of 
the spacecraft relative to a proof mass that floats freely at the center of the accelerometer and, 
second, a set of thrusters that accelerate the spacecraft. A control loop fires the thrusters 
in such a way as to null the output from the accelerometer and maintain the spacecraft on 
a purely gravitational trajectory. 



B. Gravitational Wave Response 

The response of an electromagnetic tracking signal to the passage of a gravitational wave 
has been shown by Estabrook and Wahlquist || to be a Doppler shift in the frequency of 
the received signal relative to the outgoing signal. For a gravitational wave of amplitude 
h(t), the shift in a signal of fundamental frequency v a will be given by 

= - cos 2w 

v 2 

x [(1 - cos 6)h(t) + 2 cos 9 h(t - r - r cos 6) - (1 + cos 9)h(t - 2r)] , (1) 

where r is the light travel time between end masses, 9 is the angle between the lines-of-sight 
to the probe and to the source, and ^ is a principal polarization angle of the quadrupole 
gravitational wave. 
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The angles 9 and tp will, for a general spaceborne interferometer, be slowly varying 
functions of time, dependent upon the orbital configuration of the detector. This paper 
ignores these time dependencies, which are specific to a particular mission design, so that the 
only time dependence assumed in Eq.(|l]) is the time varying amplitude of the gravitational 
wave, h(t). The sky-averaged sensitivity obtained by averaging over the antenna pattern 
associated with a particular interferometer orbital configuration will be negligibly different 
from the all-sky average obtained here for a fictitious spatially fixed interferometer^} 

It is useful to write h(t) in terms of its Fourier transform h(u). If the Doppler record is 
sampled for a time T then h(t) is related to its Fourier transform by 

,/T r +oo _ 

h(t) = — / h{u)e wt du , (2) 

where the y/T normalization factor is used to keep the power spectrum roughly independent 
of time. Using this definition of the Fourier transform, the frequency shift can be written as 

v \Pf r+°° 1 
Au(t, 9, V) = / - cos(2^) h(u, 9, 

(1 - /i) + 2fie~ iuJT(1+ ^ - (1 + fi)e- i2ujT ] e^duj , (3) 



x 



where \i = cos 9. The quantity that is actually read out by the laser interferometer tracking 
system is phase, so Eq. (§) is integrated to find the phase in cycles 

A0(t, 9, V0 = f Au(t', 9, if,) dt'. (4) 
Jo 

If the phase is divided by the laser frequency and the one-way light travel time of each arm, 
one obtains: 



2 The angles 6 and ip will vary at the orbital period of the interferometer, T or j,. In calculating the 
response to a gravitational wave of period T gw , the relative error in the power spectrum made by 
assuming the angles to be constant will be of order (T gw /T or f,) 2 . As long as the orbital period is 
much longer than the gravitational wave period, this error will be negligible. 
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x "(1 - //) + 2fxe~ iuJT{1+ ^ - (1 + fi)e' i2ujT ] —e iuJt , (5) 

where Eq. ([|) has been used to expand Au(t, 9, ■0), and arbitrary phases have been set to 
zero in the integration of Eq. (0j). The 1/u in the Fourier integral arises from the integration 
in the time domain in Eq. (01). This z(t,9,ip) has the property that, in the limit of low 
frequency, it reduces to a pure spatial strain. At high frequencies, the signal is much more 
complicated due to the three /i-terms that enter at different times. Nevertheless, in what 
follows we will loosely refer to the quantity z(t,6,ip) as the gravitational wave strain. 

The goal of gravitational wave detection is to detect the strain produced by a gravita- 
tional wave signal, as given in Eq. (^), in the presence of competing noise. The primary 
sources and spectra of this noise will be discussed in the next section. 



III. NOISE SPECTRA 
A. Common noise spectra 

The signal received in a single arm of the interferometer is given by 

Si (t) = z(t, 9 h fa) + p(t) - p(t - 2r t ) + m(t), (6) 

where Si(£) and S2(t) are the two noisy strain signals in the two arms of the interferometer, 
p{t) is the laser phase noise which is common to the two arms, rii(t) is the strain noise in 
the i th arm produced by all other noise sources, and r, is the one-way light travel time in the 
i th arm. The major source of noise in each arm, by several orders of magnitude, is the p(t) 
phase noise in the lasers. However, this noise is common to both arms of the interferometer 
and can be eliminated through signal processing. When an interferometer signal 

E(t) = si(t) - s 2 (t) 

= Zl (t) - z 2 (t) + m(t) - n 2 (t) - p(t - 2n)+p(t - 2r 2 ) (7) 
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is formed, the laser phase noise will cancel as t\ — > r 2 . Unfortunately, a space-based in- 
terferometer consisting of freely flying spacecraft will necessarily have unequal armlengths, 
preventing effective use of £(t) as a signal for data analysis. 

However, a new data reduction procedure for space interferometry has been recently 
discovered and discussed in a paper by Tinto and Armstrong |TT|] , that eliminates all laser 



phase noise even if the two light travel times t\ and r 2 are unequal. Working in the time 
domain, one defines the combination 

X(t) = Sl (t) - s 2 (t) - s x {t - 2r 2 ) + s 2 (t - 2n) 
= Zl (t) - z 2 (t) - Zl (t- 2r 2 ) +z 2 (t- 2t x ) 

+m(t) - m{t - 2r 2 ) - ri2 (t) + n 2 {t - 2n) , (8) 

which is devoid of laser phase noise for all values of the two light travel times T\ and r 2 . 
In order for the laser phase noise to cancel exactly, both light travel times must be known 
exactly so that the combination of signals in Eq. ([§]) can be correctly formed. In practice, 
however, all that is necessary is for the armlengths to be known well enough for the laser 
phase noise to become insignificant relative to the independent noise sources (discussed 
next). 



B. Independent noise spectra 

Since noise common to the two arms of an interferometer can be significantly reduced by 
the procedure just described, the limitation of sensitivity of the space-based gravitational 
wave missions is actually determined by noise sources that are independent in the two arms. 
There are two major types of such independent noise: position noise in the laser tracking 
system and acceleration noise due to parasitic forces on the proof mass of the accelerometer. 

Position readout noise is noise that acts like a change in the length of the optical path 
for the tracking signal. The strain spectral density it creates is given by 

Sn(f) = ^ 2 (9) 
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where S x is the spectral density of position noise in m 2 Hz _1 . The most important single 
source of position noise in the current designs is shot noise in the lasers. This may therefore 
be used as a rough indicator of the expected position noise for a particular interferometer. 
Shot noise produces a phase noise given by 



hv 

Pr 



(10) 



where hv Q is the photon energy and P r is the received power. The received power is calculated 
from 



Pr = Pt 



en 2 u 2 D 2 



enD 2 



(11) 



where P t is the transmitted power , the quantity in the first square brackets is the directional 
gain of the transmitter optics with diameter D and efficiency e for light of frequency z/ Q , the 
quantity in the second square brackets is the space loss at a distance r, and the quantity in 
the last square brackets is the effective cross section of the receiving optics. The effective 
position noise produced by this phase noise is 

~2 



S r 



Atx 2 u; 



;S<h. 



(12) 



Combining Eqs. (P-|T2]), the formula for the strain noise produced by shot noise in the lasers 
is 

4c 4 hu a 1 



7T 4 e 2 D A P t z/ 4 ' 



(13) 



where we have used the fact that r « cr. 

A more detailed study will include other, smaller sources of position noise. For the current 
LISA design, the strengths of many of these potential additional noise sources have been 
budgeted. The LISA shot noise is expected to have spectral density S x (f) = 1.2 x 10~ 22 
m 2 Hz -1 , while the total position noise is estimated to have spectral density S x (f) = 
1.6 x 1(T 21 m 2 Hz" 1 [§. 

In addition to position noise, there will be acceleration noise in the detectors. Accelera- 
tion noise produces a strain noise spectral density given by 
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where S a is the acceleration spectral density in m 2 s~ 4 Hz _1 . The major source of acceleration 
noise in the space detectors is expected to be parasitic forces acting on the proof mass of the 
accelerometer. These instruments are very complex and the noise is difficult to characterize, 
especially in the laboratory tests where the proof mass must be suspended in one-g rather 
than being allowed to float freely in space. Nevertheless, from what has been learned in the 
laboratory, it appears that a flat acceleration noise spectrum at a level of S a = 9 x 10~ 30 m 2 
s~ 4 Hz -1 should be achievable over most of the frequency band of interest (this is the level 
assumed in the present LISA design 0). 

The total noise curves for space gravitational wave detectors are found by adding the 
various noise spectra. Using the values for S x (f) and S a (f) given above, a total spectral 
density of strain noise in the detectors can be found by adding Eqs. (|9|) and ([TJ|) in quadra- 
ture. The resulting root spectral density of the instrument noise is plotted in Figure 0. The 
quoted LISA spectral densities have been used for the acceleration and total position noises: 
S a = 9 x 1(T 30 m 2 s~ 4 Hz" 1 and S x (f) = 1.6 x 1(T 21 n^ILr 1 . 



IV. GENERATION OF SENSITIVITY CURVES 

The sensitivity of an instrument to a gravitational wave depends on the relationship 
between the amplitude of the wave and the size of the signal that eventually appears in the 
detector. It also depends on the size of the noise in the final output signal of the instrument. 
The connection between amplitude and signal is calculated in frequency space and is called 
the transfer function. If the interferometer output £(i) has a gravitational wave contribution 
to it given by 

A(t) = Zl (t) - z 2 (t) , (15) 
then the transfer function R(u) is defined by 

S A (u) = S h (u)R(u) , (16) 
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where the gravitational wave amplitude spectral density Sh(uj) is defined by 

S h {u) = \h(u)\\ (17) 
so that the mean-square gravitational wave strain is given by 

1 poo 1 poo 

(h 2 ) = Ktfdt = — jf S h (u)dw . (18) 

Similarly, the instrumental response S^(u>) is defined such that 

W)=7^r s ^ duj > ( i9 ) 

Z7T JO 

where the brackets indicate a time average and the bar over the A in Eq. ([19]) indicates that 
it is averaged over source polarization and direction. 

In the next section, the transfer function from the gravitational wave amplitude h to 
the interferometer signal A (which is the signal part of £) is worked out. As discussed 
in Section [III A| , the preferred instrumental output would actually be X(t), defined in Eq. 
(§), since it exactly cancels the laser phase noise in the data. However, X(t) is a difficult 
quantity to work with since its transfer function will depend on the particular values of the 
interferometer arm lengths, T\ and r%. One approach to simplify the calculation is to assume 
Ti = T2 (This approach is taken in Ref. ||.) At the end of the next section, we will give the 
transfer function for X(t) with this assumption made, and will show that in this limit £(£) 
and X(t) yield the same instrumental sensitivity. 



A. Gravitational Wave Transfer Function 

Previous estimates of sensitivity || have often resorted to working in the long wavelength 
approximation, beginning from the Doppler tracking signal described by Eq. (|1|), or have 
combined results for independent single arms. Here the exact gravitational wave transfer 
function is computed, without such approximations, from the gravitational wave strain. 

Consider the geometry shown in Figure ^. The vectors Li and L2 point along the arms of 
the interferometer, and the vector k points along the propagation vector of the gravitational 
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wave. The quantities 9% measure the angular separation between the arm vectors and the 
propagation vector. The value 7 is the opening angle of the interferometer, and e is the 
inclination of the plane containing k and Li to the plane of the interferometer. 

Not shown are the principal polarization vectors of the gravitational wave, which lie in 
a plane 90 degrees away from k. The polarization angles, ipi, are measured from the point 
where the plane of the i th arm and the propagation vector intersects the plane containing 
the principal polarization vector. The angle a in Figure ||] is simply the difference of the two 
polarization angles, a = ip2 — ipi- 

The average power in the interferometer is given by 

(A 2 >=lim-/ \A\ 2 dt, (20) 

T^oo I Jo 

where A is defined by Eq. (|15D . Using the definition of z from Eq. flD this can be expanded 
to yield 

1 r°° ~ 1 

(A 2 ) = J- / du h\u)— [T\M + T 2 (u) - 2T S ( W )] , (21) 

Z7T JO (LOT) z 



where 



2it Jo ( ujr y 

Ti(u) = cos 2 (2^!) /i 2 (l + cos 2 (cg>t)) + sin 2 (c<jr) 

— 2/i 2 cos(lut) cos(a;r^i) — 2//x sin(a;r) sin(u;r/xi) , (22) 
T 2 (uj) = cos 2 (2^ 2 ) A*! (l + cos 2 (cg>t)) + sin 2 (u;r) 

— 2/i 2 cos(a;r) cos^r^) — 2/x 2 sin(u;r) sin(u;r/i2) , (23) 
T 3 (u) = cos(2^) cos(2^ 2 ) n(u) , (24) 

with jii = cos 9i, and where 

T)(u), $i, 62) = [cos(u;r) — cos(uTfjLi)] [cos(a;r) — cos(co>T/i 2 )] A*iA*2 

+ [sin(a;r) — fit sin(uTfii)] [sin(u;r) — /i 2 sin(a;r//2)] , (25) 

has been defined for convenience. The expression for the power in the detector, as given 
by Eq. ([H]), is a complicated function of frequency and of the orientation between the 
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propagation vector of the gravitational wave and the interferometer. It represents the an- 
tenna pattern for a laser interferometer gravitational wave detector. At low frequencies, 
the frequency dependence drops out and the expressions simplify greatly [13 1 , but at higher 
frequencies it remains a very complicated frequency- dependent object. 

To characterize the average sensitivity of the instrument it is customary to consider the 
isotropic power, obtained by averaging the antenna pattern over all propagation vectors and 
all polarizations: 

I i-2-k r2ir 



o Jo Jo 



(A 2 ) = — / dip I de I sinO d6 (A 2 ) . (26) 



Since the variables (Oi,ipi) and (# 2 , ^2) that figure in Eq. ( pT|) each locate the same 
propagation vector of the gravitational wave, they are not independent of one another. 
However, examination of Eqs. ( f2~2j|24|) shows that Ti depends only on the {Q\,ipi) variables, 
while T 2 depends only on (6*2, ^2) • This allows the integration of Eq. (p6|) to be performed 
very easily for the 7\ and T 2 terms without converting to a common set of angular variables. 
In each case, the averaging of the Ti and T 2 terms over the (Oi^tpi) gives precisely the same 
value: 



T 1 = T 2 = — 

8tt 

1 " 



J dipi de dd\ sinOi T\ 

A 

(27) 



1 + cos 2 (^r)) [\- j^—A + sin 2 (wr) + — i- sin(wr) cos(cjt 
2 |_ v y \S (ujtY J (lot) 6 

The average isotropic power can then be expressed as 

(A 2 ) = - / duo h\uo)-— 2 Tx - T 3 . (28) 

7T JO (iOT) z 1 J 

To complete the integration, the function T 3 (which depends on both (#i, ipi) and (6* 2 , ip 2 )) 
must finally be expressed in terms of a single set of angular variables. One may choose to 
eliminate (82,^2) i n favor of (6i,ipi) by using conventional spherical trigonometry in Figure 
[3]. For the polarization angle, the relationship is particularly simple 

1P2 = ipi + a 
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and the integration over ipi may be carried out analytically, giving, 

T3 = — - J dipi de dd\ sin 6*i T3 = — J de d9\ sin#i cos(2ct) t)(uj, 61,82) ■ (29) 

The cost of carrying out the ip\ integration is the introduction of the angle a, which can be 
related to the {e, 6>j} variables using the law of sines in Figure [3| 

sin 7 sine 

sin a = — . MO) 

sin 62 



The function T)(oj, 9i, 62) in Eq. (|25|) has terms containing /x 2 = cos 62, which must be re- 
expressed in terms of the integration variable Q\. The relationship between 61 and 9 2 is given 
by 

cos 6*2 = cos 7 cos^i + SU17 sin 6*i cose , (31) 

where 7 is the opening angle of the interferometer, and e is the inclination of the gravitational 
wave propagation vector to the interferometer. Due to the complexity of T](oj) when Eq. 
(|31"|) is substituted into Eq. (|25|), we have not been able to calculate an expression for T3 
analytically, so it will be kept as an explicit integral. 

Using the definition of R(u) from Eq. flUf) with the average isotropic power in Eq. (p8|) , 
the gravitational wave transfer function is found to be 

R(u) = 2-^—\T 1 -T 3 

[ujt) 2 L 



1 



(cor) 



cos 2 (ujt'' 



3 (cury 



4 

+ sin 2 (cur) + - — — sin(u;r) cos(a;r) 
1 ( 

— — J de d6i sin^! ^1 — 2 sin 2 cxj 17(00) 



(32) 



It is straightforward to evaluate the remaining integral using simple numerical techniques. 
The exact transfer function, including the numerical evaluation of the last integral, is shown 
in Figure |asa function of the dimensionless quantity u = ujt. 

The high frequency structure in Figure |3] dominates the shape of the transfer function at 
frequencies greater than uo ~ 1/r. At these frequencies, the armlength r of the interferometer 
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becomes comparable to the wavelength of the gravitational wave. The extrema in the transfer 
function are amplifications due to interference of the signals in the arms. The minima occur 
at frequencies 

/ = n/2r . (33) 

The appearance of these periodic amplifications is familiar from basic interferometry. 
Taking the limit of small lot in Eq. (|32|) will yield a low frequency limit for of 

4 

Sa = - sin 2 7 S h , (34) 
5 

which is in agreement with a previous result from Hellings ||[] . 



B. Sensitivity Curves 



The sensitivity curve for a gravitational wave observatory is obtained from Eq. (16) in 



the case where the spectral density of the isotropic power, is equal to the spectral density 
of noise in the detector. The noise in the final instrumental signal may be related to the 
hardware noise in each detector by inspection of Eq. (|7|). Assuming that n\(t) and ^(i) are 
uncorrelated, the spectral density of the signal produced by phase measurement noise is 

Sn = 4S n , (35) 

where S n is the spectral density of noise in a one-way measurement. One factor of 2 in Eq. 
(|35|) comes from the two rii(t) in Eq. (|7|) adding in quadrature, since they are uncorrelated, 
and one factor of 2 arises from the fact that the two-way noise in each arm is a/2 times the 
one-way noise. The final equation for gravitational wave amplitude sensitivity is thus 



3 In the notation of this paper, the interferometer signal has been multiplied by the laser period 
and divided by the arm length of the interferometer, as described in Eq. (||). To agree with the 
expressions for the spectral density of the interferometer signal Sg in the literature, this factor 
must be accounted for above, so that S$ = {tv ) 2 S^. 
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Sh = li = li = ir- (36) 
The spectral amplitude sensitivity is simply the square root of Sh, or 

h S = 4K = 2^. (37) 

Using the noise curve given in Figure |2|, and using the LISA value of cr = 5 x 10 9 m, the 
LISA sensitivity curve, computed using Eq. (0), is shown in Figure ||. The position of the 
high frequency "knee" occurs at / = 1/(2ttt) = 10~ 2 Hz. 

The greatest sensitivity is seen to occur in a mid-frequency "floor" ; the level of this floor 
is set by the size of the position noise. The width of the floor is a function of the acceleration 
noise level and the arm length of the interferometer. The low frequency rise occurs when 
acceleration noise begins to dominate over position noise; the high frequency rise is caused 
by the turnover in the transfer function, at / ~ l/(27rr). 

C. Comparison with Prior Results 



One may compare the transfer function of Eq. (|3~2j ) with others that have previously 
appeared in the literature. These transfer functions differ from the one presented in this 
work in two ways. 

First, previous results in the literature have approximated the transfer function either by 
working in the low frequency limit, where the transfer function becomes constant or by 
combining single arm results rather than treating the interferometer as a single instrument 
0. In the latter cases, the transfer functions do not include the interference term between 
the two arms of the interferometer (the term denoted T 3 in Section [IV A|) . 



4 The full transfer function has not been needed for the ground-based interferometers such as 
LIGO since in their operation (as Fabry-Perot cavities) the low frequency limit is valid up to about 
40 - 50 kHz. 
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Second, the transfer functions are often incorrectly normalized [0,0. The problem is 
that the transfer function between gravitational wave amplitude and detector response is a 
well-defined quantity (see Eq. (|IED), which is inappropriate to renormalize. When averaged 
over position and polarization, the value of R(u) is given by Eq. fl32|). In the low frequency 
limit this becomes R{uj) = (4/5) sin 2 7 = 3/5 (Eq. fl3~3| ) with 7 = 60°). However, the transfer 
function as given in Refs. |||7| has been normalized by dividing by the maximum sensitivity 
of the instrument, i.e., by the response for a source with optimal polarization located at a 
point normal to the plane of the detector. As may be seen by substituting /ij = in Eqs. 
( p| - |25|) , this maximum instrument response is (A 2 ) = 3(h 2 ), so division by this quantity 



would give a transfer function whose low- frequency limit is R = 1/5 rather than R = 3/5. 
The previous LISA sensitivity curves thus underestimate the strength of the signal that a 
gravitational wave of amplitude h will produce in the detector (see Eq. (16)) and conclude 
that the minimum detectable amplitude is 

h L i S A = V3h , (38) 

thereby underestimating the sensitivity of the detector. 



D. The Transfer Function for X(t) 

As discussed at the beginning of Section |IV|, the preferred signal from the interferometer 
for purposes of data analysis is not but X(t), since the common laser phase noise exactly 
cancels in X(t). Unfortunately, the transfer function for X(t) depends in a complicated way 
on the particular values of 77 and r 2 . However, substantial simplification occurs in the special 
case 7~i ~ T2 = r, which was treated in Ref. |J. In this case, the formula for X(t) becomes 

X(t) = Zl (t) - z 2 (t) - [zi(t - 2r) - z 2 (t - 2t)] 

+m(t) - n 2 (t) - [m(t - 2r) - n 2 (t - 2r)] , (39) 

which may also be written as 
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X(t) = E(t) + a(t) (40) 

where 

S(t) = A(t) - A(t - 2r) , (41) 

and 

a{t) = N(t) - N(t - 2r) . (42) 

Here A(t) is given by Eq. (|T§) and N(t) = m(t) - n 2 {t). 

The portion of X(t) which is a gravitational wave signal is then, by Eq. fl41~l) , simply two 
copies of the A signal we have previously analyzed. The transfer function for H(i) may be 
calculated following the procedure in Sec. |[TV Aj| , to find 



Ss(uj) = 4 sin 2 (cot) S^(uj) 

= 4sm 2 (uoT)S h (uj)R(uj) , (43) 

where Eq.([R]) has been used for the second equality. The transfer function for X will then 
be 

Rx(uj) = Asin 2 (ujT)R(uj) , (44) 



where R(u) is the previously derived transfer function for A, given by Eq. (|32D. 

At first glance, it appears the that transfer function for X could be as large as four times 
the transfer function for A, and hence that weaker gravitational waves could be observed 
by constructing X. However, to determine the sensitivity of the interferometer, one must 
equate Ss(a>) in Eq. (fHf ) to the appropriate spectral density of noise. For X, this is formed 
using the a(t) combination defined in Eq. (|42|). Because the noise contribution to X is 
formed by the same subtraction process, it will similarly have 

S a (oo) =4sin 2 (ujT)S N (uj) . (45) 

The additional factors due to the subtraction process used to form X, which appear in Eqs. 



45|), will cancel in the computation of the sensitivity limit of the interferometer: 



_ 4 sin 2 (lot) S n (u) _ 4S^ ^ 
Rx R 

where the second equality follows from Eq. (|44]) and the definition of N in terms of the n^. 

Thus, in the limit as n — > r 2 , the sensitivity curve for X(t) will be identical to that 
previously computed for A(t), as shown in Figure |5|. We are investigating transfer functions 
for unequal arm cases and intend to present these in a future paper. 

V. DISCUSSION AND RECONCILIATION 

The final goal of producing noise curves and response functions is to answer the question 
of what gravitational wave sources are detectable. Several different approaches have been 
used in the literature to answer this question. It is the purpose of this section to discuss 
these different approaches and to try to reconcile them. 

Section |T| discussed the actual strain noise spectral density (S n (f), see Figure 0), and 



Section |IV] showed how the gravitational wave transfer function (R(f), Figure |j) is used to 
determine the gravitational wave sensitivity curve (Sh{f), Figure |5|). The gravitational wave 
transfer function actually represents an averaged response of the interferometer, averaged 
over waves coming from different directions and having different wave polarizations. The 
transfer function, however, ignores the nature of the source creating the wave; no averaging 
was performed over parameters describing the sources themselves. Thus, the resulting grav- 
itational wave sensitivity curve is equally appropriate for most types of sources. Of course, 
because of the averaging over direction, it is not completely appropriate for a source (like 
the interacting white dwarf binary AM CVn) whose direction is known. Nevertheless, this 
sensitivity curve seems to us to be a valuable intermediate tool to characterize the capability 
of gravitational wave detectors, especially since the strengths of the waves from many differ- 
ent types of sources may be plotted on this same graph. To calculate the gravitational- wave 
strength for a particular type of source as it is to be plotted on the sensitivity curve graph 
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requires some additional steps specific to each class of sources. 



A. Sensitivity to Particular Types of Sources 

For continuous monochromatic sources like circular compact binaries, the ultimate source 
sensitivity comes when the source is sampled for a long period of time so as to narrow the 
bandwidth, Af. A continuous source with frequency / and amplitude h that is observed 
over a time T will appear in a Fourier spectrum of the data as a single spectral line with 
root spectral density 

h f = = hVf. (47) 

The frequency / and the amplitude h will depend on the several parameters of the binary 
such as total mass, semi-major axis, distance to the binary, inclination to the line-of-sight, 
etc. If one wants to consider a population of monochromatic binaries, then an average 
over the parameters of possible systems should be taken, and the value used for h at each 
frequency should reflect this process. 

For signals that are short bursts, such as would be produced by a compact body in highly 
elliptical orbit about a massive black hole or by high-velocity encounters of massive compact 
bodies in a region of high stellar density, the signal would typically have a characteristic 
pulse width r with a broad spectral content and would occur only once in the lifetime of the 
detector. The detection of such a source is optimal when the bandwidth Af is not much 
larger than r _1 . The relationship between the amplitude of the pulse h and the strength of 
the signal as it would be plotted on the root spectral density graph is therefore given by 

hf = ^j = hV^- (48) 

Again, if a population of burst sources is to be considered, then one should average the 
amplitudes over the possible parameters of the source, including all source orientations and 
accounting for the short duty cycles that such sources typically have. 
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Finally, one may consider how to include a stochastic background of gravitational waves 
as a source on the gravitational wave sensitivity graph. In some cases, such as the case of a 
stochastic superposition of close compact binary stars, the amplitude spectral density Sh is 
known. In these cases, the root spectral density 

hf(f) = y/sitf) (49) 

can be plotted directly on the sensitivity graph. In other cases, what is assumed is a spectrum 
of the cosmic energy density, Q(f). The relationship between amplitude and energy density 
is given by 

s h (f) = ^n(/) (so) 

where H is the Hubble constant. The square root of Sh{f) may again be plotted directly on 
the sensitivity graph. Finally, a cosmic background is often assumed to be simply peaked 
at some frequency, f p , and spread over a decade of bandwidth. In this case, the bandwidth 
will be approximately equal to the peak frequency, and a point plotted at 

1^1 2 Jp' 

would represent the peak of the decade-wide spectrum corresponding to a total integrated 
energy density of Q. 

B. Reconciliation with Other Results from the Literature 

Substantially different methods have often been used in the previous literature to il- 
lustrate the capability of a spaceborne interferometer to detect astrophysical sources of 
gravitational waves. This section describes how to relate the methods used in this work to 
the approach usually taken in the literature associated with the proposed LISA mission || , 
and the comprehensive development of sensitivity limits provided in the review article of 
Thorne (T7| . 
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As in this paper, other analyses create an initial sensitivity curve for LISA by analyzing 
the various noise sources in the interferometer and then dividing by a transfer function [[|[7| . 
The LISA instrumental sensitivity, in terms of hf, is generally converted to an effective 
gravitational wave amplitude sensitivity, h, by dividing the value of hf at each frequency by 
the square root of an assumed one year integration time. The resulting sensitivity curve in 
h is also multiplied by a factor of 5, so that the final curve indicates a one-year integrated 
threshold for a signal with signal-to-noise ratio > 5. The combination of these factors is 

/4 y NR 5 = 8-9 x 10~% , (52) 

where hf is the spectral amplitude in units of Hz -1 / 2 . This relation assumes that the spurious 
factor of \/3 noted in Section [IV Q has been corrected. 

It should be emphasized that, after having assumed a one-year integration time, it is only 
appropriate to plot monochromatic gravitational wave sources against this type of sensitivity 
curve, and only for a case where a signal-to-noise ratio of 5 is actually required. The danger 
in using such a graph to characterize the overall detector sensitivity is that it is subject to 
easy misinterpretation, which has often occurred in the literature. Examples include Hogan 
14j , which overestimates the LISA sensitivity to a cosmic stochastic background by over 



three orders of magnitude by misunderstanding how the 1-year integrated curve would apply 



to a stochastic signal, and Aguiar et. al. [15|, which significantly overestimates the LISA 



sensitivity to bursts from black hole oscillations lasting a few minutes by using the curves 
representing coherent integration of a signal for a year.Q On the other hand, the usual LISA 
sensitivity curve graph underestimates LISA's sensitivity to monochromatic known sources, 
such as AM CVn, for which a signal-to-noise ratio of 5 is an excessive requirement. 

Another approach is taken in the comprehensive review article by Thorne |T7|]. In this 



review, Thorne carefully distinguishes between the three types of source, and analyzes the 



5 The LISA team has recognized the difficulty of using such plots to characterize the detector 
sensitivity to bursts [16]. 
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response of a detector to each type. For each type of source, he finally derives a detector 
sensitivity which is denoted by h 3 / yr (f), defined as: 

h/yr(f) = 11 [fS h (f)] 1/2 burst sources, (53) 
h/yr(f) = 3.B[S h (f) x l(T 7 Hz] 1/2 periodic sources, (54) 

h/yAf) = 4-5 f 1Q _7 Hz 1 [fSh(f)} 1/2 stochastic background . (55) 

Here Sh(f) is the same spectral density of detector noise utilized in this paper. The derived 
sensitivity, /i3/ yr , represents the weakest level at which one has 90% confidence of detecting 
a signal using two cross-correlated identical detectors. For burst sources, the level is set by 
demanding three detections per year. For periodic sources, h 3 / yr represents the level obtained 
in a 1/3 year integration, assuming the frequency and phase of the source are known. Finally, 
for stochastic sources, h 3 / yI represents the weakest source that can be detected with a 1/3 
year integration of the cross-correlation function between two independent detectors. 

The development of Thome's three dicta for detection of burst, periodic, and stochastic 



sources involves detailed assumptions which are described at length in Ref. |[L7 |; that dis- 
cussion will not be reproduced in full here. It is worth noting, however, that not all of the 
assumptions made in the development of these dicta are valid for a space-based interfer- 
ometer. For example, Thorne assumes one has two independent cross-correlated detectors. 
While this is a valid assumption for LIGO, it is not for a LISA-style instrument. Although 
the data from the three LISA arms can be treated as two interferometers, and hence mea- 
sure both polarizations of gravitational waves, the two interferometers are not completely 
independent, as they share an arm. 



VI. SUMMARY 



Constructing the sensitivity curves for a spaceborne observatory, such as those shown in 
Figure ||, is the first step in understanding the response of the instrument to gravitational 
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radiation. The formalism developed here can be used to determine the sensitivity of any 
space-based interferometer simply in terms of the essential parameters which describe the 
overall design of the instrument. Not only does this allow one to quickly and accurately assess 
the performance of one proposed observatory compared to another, but it also provides a 
quick and easy method for considering new observatory designs. 

Carefully detailing the response of any new instrument prepares us for the inevitable 
detection of unexplainable signals from distant astrophysical sources, and provides a clear 
idea of how to improve our instrumentation for the construction of the next generation of 
observatories. 
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FIGURES 

FIG. 1. A typical configuration for a spaceborne gravitational wave interferometer. Three 
probes form an equilateral triangle, inscribed on the relative orbits of the spacecraft. A 
single Michelson interferometer is formed using any two legs of the triangle, and a second 
(non-independent) interferometer can be formed using the third arm of the configuration in con- 
junction with an arm already in use for the primary signal. 

FIG. 2. The root spectral density of the noise in the LISA interferometer formed by adding the 
strain noises induced by acceleration noise and by total position noise in quadrature. 

FIG. 3. The geometrical relationship of the interferometer to the propagation vector of a 
gravitational wave, used to conduct spatial averaging. The arms are designated by vectors Li and 
L2 (solid black vectors), while the propagation vector of the gravitational wave is given by k (open 
white arrow). One arm of the interferometer is aligned along the polar axis of a 2-sphere, the other 
arm lying an angular distance 7 away long a line of constant longitude. The angles 6{ relate the 
vector k to the arms of the interferometer, and the angle e is the inclination of the plane containing 
k and Li to the plane of the interferometer. 

FIG. 4. The transfer function R{u) is shown as a function of the dimensionless variable u = lot. 
Note that it is roughly constant at low frequencies, and has a "knee" located at u = lot ~ 1. 

FIG. 5. The sensitivity curves for the proposed LISA observatory is shown. The low frequency 
rise is due to acceleration noise in each of the systems. The high frequency rise is due to the "knee" 
in the transfer function at / ~ (2ttt) . The structure at high frequencies is a consequence of the 
high frequency structure in the gravitational wave transfer function. 
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